Measuring spin and charge correlations via tunneling-current conductance fluctuations 
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Scanning tunneling miscoscopy is one of the most powerful spectroscopic tools for single-electron excita- 
tions. We show that the conductance fluctuations, or noise in the conductance, of a tunneling current into an 
interacting electron system is dominated by density-density and spin-spin correlations. This allows one to probe 
two-particle properties (susceptibilities) and collective excitations by standard experimental tunneling methods. 
We demonstrate this theoretically, using a novel many-body calculation for the multi-center Kondo problem, in- 
cluding both direct and indirect exchange between magnetic atoms. An example of the two-particle correlations 
around a single magnetic adatom in the Kondo regime, as would be viewed by a scanning tunneling microscope, 
is given. The spatial dependance of the charge and spin correlations, including the formation of the Kondo cloud 
in the spin sector, are shown. 
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Noise spectroscopy, typically current noise, has become an 
exceedingly useful tool in the study of electronic systems. The 
intrinsic noise, i.e. the fluctuations of a signal due to inher- 
ent uncertainties, generated by an electronic system is not a 
simple set of random uncorrelated events but contains funda- 
mental information about electron-electron correlations, that 
is not seen in the (averaged) signal itself. Perhaps, the most 
well-known demonstration of this was the experimental ver- 
ification of fractionally charged quasi-particles in the quan- 
tum Hall regime, by shot-noise measurements 1 ' 2 . The study 
of current-current correlations, such as shot-noise or Johnson- 
Nyquist (thermal) noise, has a long history and has been a 
prevalent topic of both experimental and theoretical investi- 
gations, along with work in related areas, such as conduc- 
tance fluctuations of mesoscopic wires, mostly by studying 
the effects of disorder 3 ' 4 ' 5 . Here, we show that in the weak- 
local-tunneling limit the conductance fluctuations of an unpo- 
larized (spin) current into an interacting system is dominated 
by density-density correlations or for a spin-polarized current 
by the spin-spin correlations. This allows one to extract two- 
particle characteristics from tunneling experiments along with 
single-particle quantities, such as the density of states or mag- 
netization. Macroscopic two-particle properties, like the com- 
pressibility or magnetic susceptibility, are easily measurable. 
However, few if any techniques exists to measure these quan- 
tities at a local microscopic scale, with spatial and energy res- 
olution. For instance, almost fifty-years after the discovery 
and theoretical explanation of the Kondo effect, the spatially 
localized spin correlations around a magnetic impurity — the 
Kondo cloud — has never been experimentally observed. 

The atomic spatial resolution of a scanning tunneling mi- 
croscope (STM) makes it a natural choice to study local fluc- 
tuations on the microscopic scale. This combination of an 
STM and noise spectroscopy has already been used to de- 
velop the field of electron spin resonance scanning tunneling 
microscopy (ESR-STM) 6 . With a similar experimental setup 
in mind we will limit ourselves to an STM system, along with 
making the common approximations, such as weak local tun- 
neling and where the STM is assumed to be a weakly or non- 
correlated Fermi liquid with a featureless single-particle den- 



sity of states. Although, the work presented here is with ref- 
erence to an STM, the results are not limited to such systems 
and are valid anywhere such approximations can be made. 
The total Hamiltonian is taken to be 

H = H STM + H suh + H tnn , 

where H suh is the general interacting Hamiltonian of the sub- 
strate and H STM = £ k(T (e k(T - (i- eV)al a a krJ is the Hamil- 
tonian of the STM. As usual, we assume that the STM is 
a noninteracting Fermi system with an energy independent 
density of states near the Fermi energy e F = /i. This is an 
optimal condition for spectroscopic aims. The chemical po- 
tential of the STM is displaced by eV, where the charge of 
the electron is — e and V is the applied voltage. The tunnel- 
ing is determined by H tun = £k,k',a [Ta k Av + Hc ]> 
with tunneling amplitude T, which is assumed to be inde- 
pendent of momenta k,k'. The operators b^, a and b k , a are 
the mode creation and annihilation operators for the substrate. 
The current operator is defined as I = —ed t NsTM, where 
^Vstm = J2ka a ko- a k(T i s tne particle number operator for the 
STM. Assuming [Astm, H suh \ = 0, by Heisenberg's equa- 
tion of motion, I = ie[iVsTM, -^tun] (with ft = 1). Evaluat- 
ing the commutator using {a k(T , 6 kv , } = leads to, within 
the tunneling Hamiltonian formalism, the common expression 
for the current operator 7 

^ K EE[ Ta Uv-H.c.]. (i) 

a k,k' 

To obtain the experimentally measured current, one needs to 
obtain the non-equilibrium expectation value of (1). We do 
this within the linear response (LR) regime, treating the tun- 
neling as the perturbation and assuming the STM and sub- 
strate are separately in thermodynamic equilibrium. If the 
system is decoupled in the infinite past, t = — oo, the current 
operator within LR is given as 

t 

I hR (t) = i + i J dt'[H tun (t')J(t)], (2) 



where 6{t) = e iH °*6 e - iff °* and H = H STM + H suh . 
The expectation value of (2) with respect to Ho, (l LR ) H = 

Tr I LR e~P H ° /Tr e~ l3H o, gives the current to leading order 
in the tunneling amplitude, T. Therefore, the linear conduc- 
tance, in this approximation, is given by G LR = d v (l LR ) H . 
Assuming that the expectation value and derivative commute, 
a conductance operator can also be defined by G LR = d v I LR 
(see supplemental material for details). With an operator ex- 
pression for the conductance, one can obtain the fluctuations 
or specifically the spectral density of the conductance, defined 
as 

S(r,u) = ^ J dt e iu)t ( { <5Glr (r , t) , 5 Glr ( v i^)})h ' ^ 

where <5Glr = Glr — (Glr). The low-temperture zero- 
frequency limit of (3) can be shown to be given by 

S{t,w = 0) = ^ e 4 |T| 4 [p STM ( e l/)] 2 x^ b (r,c - 0) 

+ 327r 2 e 4 |T| 4 [m STM (eV)] 2 ^ b (r, u = 0), 

(4) 

With psTM = E CT PSTM' ™STM = PsTM ~ PsTM' and 

X £ h (r,t) = (5h(r,t)5n(r,0)) Hgah (5a) 
xSb(r,t) = (s z (r,t)F(r,0)> ifsub , (5b) 

where Xsub and Xsub sire trie l° ca l charge- and spin- 
susceptibilities respectively, with density and spin-density op- 
erators of the substrate n and s z . The compressibility (5a) is 
given in terms of the density variation; Sn = n — (n). The 
extension of equation (4) to non-zero frequency is in principle 
straightforward, although relating the finite-frequency results 
to physically meaningful quantities is not. This is analogous 
to the standard current shot-noise result, where it is only the 
zero-frequency component of the noise that is directly pro- 
portional to the measured current. Equation (4) in itself may 
not be surprising as, loosely speaking, the conductance is de- 
termined by a single-particle correlation function (the density 
of states), thus for a conductance-conductance correlation one 
could expect two-particle quantities. As a result one can ob- 
tain local susceptibilities, as a function of position, using only 
a single STM or current probe. 

It should also be noted that, in the weak tunneling limit 
the standard expression for the current noise, which is given 
by the current-current correlation function, leads to the well- 
known shot-noise relation. The zero-frequency shot noise is 
proportional to the current itself which goes as the tunneling 
amplitude squared, while (4) is proportional to |T| 4 . This 
seemingly contridictory result is explained by the fact that in 
general the complete characterization of fluctuations or noise 
of any signal is not determined solely by a second-order mo- 
ment, such as a current-current correlation (variance), but by 
all higher moments as well. For instance one could obtain 
a similar result to equation (4) from the current signal itself, 
by a suitable choice of a current-current-current-current cor- 
relation (kurtosis). Also, in any real experimental measure- 
ment, the tunneling amplitude enters the tunneling current 
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FIG. 1 : STM setup. Schematic of an STM and several magnetic im- 
purities with spin S on a non-magnetic surface. Exchange between 
the magnetic atoms can be indirect, though the substrate, or if in close 
proximity, directly by the overlap of impurity wave functions. Nor- 
mally, the Kondo effect is probed by tunneling into a single impurity 
or nano-cluster, i.e. r w 0, where the formation of a Fano-lineshape 
in the current(conductance)-voltage curve is observed. 

to all orders. These higher order terms, the so-called ver- 
tex corrections, to the tunneling are small, and they are typ- 
ically neglected; although recent STM experiments 8 ' 910 have 
shown they can lead to detectable contributions. For exam- 
ple, the single-triplet transition (which appears in the spin- 
susceptibility but not in the single-particle density of states) 
of an atomic spin chain has been observed. As yet no full the- 
oretical description of these effects exists, for one has to go 
beyond the standard tunneling Hamiltonian approximations 1 1 
to describe the experimental dy(l(eV)) curves. Such a for- 
mulation, relating these vertex corrections to physical quan- 
tities, like spin or charge susceptibilities, is highly desirable 
and is of ongoing interest. 

As can be seen from (4), for a non-magnetic STM, i.e 
jbstm = 0, the charge susceptibility (5a) determines the fluc- 
tuations, while for a spin-polarized STM (SP-STM), tostm 7^ 
0, the spin susceptibility (5b) would be expected to domi- 
nate. This is analogous to standard SP-STM measurements, 
where the magnetic structure of the substrate is easily re- 
solved, even for relatively small spin polarization 1 2 ' 1 3 . Among 
the many possibilities, one could spatially resolve the low- 
energy spin correlations near one or more Kondo impurities, 
where the geometry of the nano-cluster, as well as the direct 
and indirect exchange processes between atoms are in strong 
competition 14 with the Kondo correlations, see Fig. 1. Here, 
as a simple but intriguing example, we show that one can ex- 
plore the localized spin correlations around a single magnetic 
atom, the so-called Kondo cloud. We now turn to a calculation 
of the local susceptibilities (5a) and (5b) for such a system. 

The Kondo effect has been extensively studied both theoret- 
ically and experimentally 15 , more recently by using an STM 
to image single or multiple magnetic adatoms on a metallic 
surface 16 ' 17 ' 18 ' 19 ' 20 ' 21 , e.g. Fig.l. Experimentally, for the most 
part, attention has been restricted to measuring the forma- 
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tion of the Abrikosov-Suhl-Kondo resonance of the density of 
states, while in the Kondo regime. Theoretically many other 
quantities have been explored, such at the non-local correla- 
tions between the impurity and conduction electron, which is 
typically used to define and determine the Kondo screening 
length £ K . To calculate the required bath correlation func- 
tions (5a) and (5b), in the presence of magnetic impurities, we 
used a numerically exact scheme, briefly outlined below (see 
supplementary information for details). In principle this for- 
malism allows one to calculate all physical quantities includ- 
ing those of the impurities, the bath, and all n-particle corre- 
lations. Although straightforward and when applied recovers 
known results 22 ' 23 ' 24 ' 25 , to our knowledge it has not appeared 
in the literature. We believe this is the ideal method to ap- 
ply to real experimental systems. Especially those involving 
multi-impurities, in close proximity, where exchange between 
atoms, direct and indirect, or even the overlapping and inter- 
ference of individual Kondo clouds, becomes important. 

The most general Hamiltonian of a noninteracting substrate 
with n-bands coupled to A-atomic impurities with amplitude 
^ko- as' anc ' including direct exchange J a , a > between impuri- 
ties is 

N 

H suh = ^ J2( £ Wa ~ M)&n,kA,k CT + Y Y E aS ci s C as 



n k,cr 



a—1 s 



2 y , y t ^s 1 ---s i c as 1 c as 2 C as 3 c as 4 + r> ^ t Ja.a'^a ■ S a 
a—1 S1---S4 

N 



(6) 



a — 1 n k.cr,s 



where c' as (c as ) is the electron creation (annihilation) oper- 
ator for an impurity, with a complete set of quantum num- 
bers s. Here S a is the total spin of an adatom, E as are the 
bare energy levels, and Uf the Coulomb interaction. In 
principle all of the above parameters, along with the disper- 
sion of the metal e^ CT , could be obtained from an ab initio 
calculation, e.g. density functional theory. With respect to 
the Hamiltonian, equation (6), the generating functional, with 
action Sa, for the entire system can be written as a func- 
tional integral over Grassmann variables, including source 
terms Ai and Ai for the bath electrons and each impurity; 
Z4 = J D[c, c]D\b, b]e~ SA . Because the host metal is as- 
sumed to be noninteracting, i.e. Gaussian, the bath electrons 
can be integrated out exactly, leading to a reduced generat- 
ing functional Za ~ f D[c,c]e~ SA , with an effective ac- 
tion (Hamiltonian) for the impurity sites. The propagator of 
the bath or any correlation function can be obtain by suitable 
functional differentiation of the effective action with respect 
to the sources. Doing so, the only unknown correlators are 
those of the impurities. The evaluation of which can be done 
using a variety of computationally fast and accurate impurity 
solvers. Here, we used the numerically exact continuous-time 
quantum Monte Carlo (CT-QMC) method of Ref. [26]. 

As an example of the usefulness and flexibility of the 
above formalism, we calculated the zero-frequency non-local 
charge and spin correlations between an impurity and the bath 
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FIG. 2: Non-local impurity-bath correlations. The zero- 
frequency non-local (NL) spin (top) and charge (bottom) correla- 
turns, denned as XN P L ( r .*) = <s z (r, t)5f mp (0)) and X N L (r, t) = 
(<5n(r, t)8Ni mp (0)\, between a magnetic impurity and the conduc- 
tion electrons as a function of distance from the impurity. In the 

2 1 k I ^ I r I ^ 

graph xnl — Viv r ^ NL • AH other relevant parameters are given 
in Table I. 
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TABLE I: Parameters for the symmetric single-impurity spin-1/2 
Anderson model in energy units of F — w\V\ 2 N a (eF), where ef 
and JV<T(ep) are the Fermi energy and spin-resolved noninteracting 
density of states of the bath respectively. Lengths are in units of the 
inverse Fermi wavevector kp | ~ 1 . The Kondo temperature Tk is ob- 
tained from the Bethe-Ansatz solution" 7 . The expected size of the 
Kondo screening cloud is of order £ K m hvF/(k B T K ), where vf 
is the Fermi velocity. Unless otherwise stated all calculations were 
done at an inverse-temperature of (fc B T) _1 = j3 = 200, well below 
the Kondo temperature for these parameters. 



(Fig. 2), which has been extensively studied both analytically 
and numerically for equal times 23 ' 24 ' 28 ' 29 ' 30 , i.e. t = 0. For 
simplicity and clarity, the real Hamiltonian (6) is approxi- 
mated by the symmetric single-impurity spin-1/2 Anderson 
model 15 , with an onsite U and a 3D parabolic dispersion for 
the bath. We will also neglect the direct tunneling into the im- 
purity. With this contribution our results would be modified 
only for r w 0. In the low-energy or low-frequency regime, 
Fig. 2 clearly shows a separation of scales between spin and 
charge correlations as the on-site Coulomb energy U is in- 
creased, as compared to the case of equal times (effectively 
high-energy), e.g. Ref. [28]. This is what one would expect as 
the dominate correlations of the Kondo effect are at low en- 
ergy, while charge fluctuations of the impurity are suppressed 
with increasing U. 

In Fig. 3 the zero-frequency local charge and spin suscepti- 



4 



4 



h 1 



1 1 1 1 1 1 1 1 1 1 — 




1 1 1 1 1 1 1 1 


U = 1 




U = 2 











5 10 15 20 25 



| 1 1 r | 

FIG. 3: Local correlations of the bath. The zero-frequency lo- 
cal spin (top) and charge (bottom) susceptibilities, equations (5 a) 
and (5b), of the bath for the single-impurity spin-1/2 Anderson 
model as a function of the distance from the impurity site. Here, 
X = [ 2r |'^^ r ^ ] X- All parameters are the same as in Fig. 2 and are 
given in Table I. 




FIG. 4: Envelope of spin-spin correlation. The envelope of the 
zero-frequency local spin susceptibility (equation (5b), shown in 
Fig. 2)for fixed U = 1 and different temperatures T; (3 = (fc B T) _1 . 
The horizontal axis has been rescaled by the Kondo length; £ K = 
Hvf I (fc B T K ). All other parameters are given in Table I. 



bilities of the bath [equations (5a) and (5b)], which determine 
the conduction fluctuations, for different U are shown. The 
oscillation period of both is ir/2, the same as Friedel oscil- 
lations, but the spin correlation function is ir/2 phase shifted 
from that of the charge. The envelope or decay of the spin 
correlations as a function of distance is shown in Fig. 4. For 
r <C Ck tne correlations show a non-algebraic decay, which 
changes to a power-law for r sa £ K , and ultimately, at fi- 
nite temperature, the correlations are exponentially cut-off 
by the thermal length £ T ~ hv F /k B T. The appearance of 
the thermal length can be seen for the largest values of r in 
Fig. 4, as the power-law changes over into an exponential. 
At zero-temperature the decay would remain a power-law for 
r 3> £ K . The simple physical interpretation of this is: at zero- 
temperature the magnetic impurity is almost fully screened 
by conduction electrons within £ K , thus outside of this length 
scale correlations are weak and decay quite rapidly, but within 
£ K correlations of the bath, mediated by the impurity, remain 
non-trivial. 

In conclusion we have shown that the conductance fluctu- 
ations of a tunneling current into an interacting system is de- 
termined by the charge and spin susceptibilities of the sys- 
tem. We have also shown that one application of this is to use 
an SP-STM to detect the Kondo screening length, £ K , for a 
magnetic adatom on a metallic surface. Although the Kondo 
problem has been and continues to be one of the most inten- 
sively studied phenomena in condensed matter physics, there 
has yet to be an experimental conformation of this theoretical 
prediction concerning the screening cloud. We have further- 
more developed a general method to exactly calculate n-point 
correlations for experimentally relevant setups consisting of 
multiple adatoms or correlated "sites". Extension of these re- 
sults to superconducting systems and quantum dot geometries 
is of future interest. 



Acknowledgments 

This work has been supported by the German Research 
Council (DFG) under SFB 668. K.P. would like thank Kirsten 
von Bergmann, Germar Hoffmann, Vaclav Janis, Hartmut 
Monien, and Jens Wiebe for useful discussions. 



* kpatton@physnet.uni-hamburg.de 

1 Kane, C. L. and P. A. Fisher, M. Nonequilibrium noise and frac- 
tional charge in the quantum Hall effect. Phys. Rev. Lett. 72, 724 
(1994). 

2 Saminadayar, L., Glattli, D. C., Jin, Y. and Etienne, B. Observa- 
tion of the e/3 fractionally charged Laughlin quasiparticle. Phys. 
Rev. Lett. 79, 2526 (1997). 

3 Lee, P. A. and Stone, A. D. Universal conductance fluctuations in 
metals. Phys. Rev. Lett. 55, 1622 (1985). 

4 van-Oudenaarden, A., Devoret, M. H., Visscher, E. H. , Nazarov, 
Yu. V. and Mooij, J.E. Conductance flucuations in a metallic wire 
interrupted by a tunneling junction. Phys. Rev. Lett. 78, 3539 
(1997). 



5 von-Oppen, F. and Stern, A. Electron-electron interaction, con- 
ductance fluctuations and current noise. Phys. Rev. Lett. 79, 1114 
(1997). 

6 Manassen, Y, Hamers, R. J., Demuth, J. E. and Castellano, A. 
J. Direct observation fo the precession of individual paramag- 
netic spins on oxidized silicon surfaces. Phys. Rev. Lett. 62, 2531 
(1989) 

7 Mahan, G. Many Particle Physics (Plenum Publishers, New York, 
2000 ) 

8 Hirjibehedin, C. F, Lutz, C. P. and Heinrich, A. J. Spin coupling 
in engineered atomic structures. Science 312, 1021 (2006). 

9 Heinrich, A. J., Gupta, J. A., Lutz, C. P. and Eigler, D. M. Single- 
atom spin-flip spectroscopy. Science 306, 466 (2004). 



5 



10 Heinrich, A. J., Lutz, C. P., Gupta, J. A. and Eigler, D. M. 
Molecule cascades. Science 298, 1381 (2002). 

1 1 Prange, R. E. Tunneling from a many-particle point of view. Phys. 
Rev. 131, 1083 (1963). 

12 Wortmann, D., Heinze, S., Kurz, Ph., Bihlmayer, G. and Bliigel, S. 
Resolving complex atomic-scale spin structures by spin-polarized 
scanning tunneling microscopy. Phys. Rev. Lett. 86, 4132 (2001). 

13 Heinze, S., Bode, M., Kubetzka, A., Pietzsch, O., Nie, X., 
Bliigel, S. and Wiesendanger, R. Real-space imaging of two- 
deimensional antiferromagnetism on the atomic scale. Science 
288, 1805 (2000). 

14 Savkin, V. V., Rubtsov, A. N., Katsnelson, M. I. and Lichtenstein, 
A. I. Correlated adatom trimer on a metal surface: a continuous- 
time quantum Monte Carlo study. Phys. Rev. Lett. 94, 026402 
(2005). 

Hewson, A. C. The Kondo Problem to Heavy Fermions (Cam- 
bridge University Press, Cambridge England, 1993). 

16 Li, J., Schneider, W., Berndt, R. and Delley, B. Kondo scattering 
observed at a single magnetic impurity. Phys. Rev. Lett. 80, 2893 
(1998). 

17 Madhavan, V., Chen, W., Jamneala, T., Crommie, M. F. and 
Wingreen, N. S. Tunneling into a single magnetic atom: spectro- 
scopic evidence of the Kondo resonance. Science 24, 567 (1998). 

18 Knorr, N., Schneider, M. A., Diekhoner, L., Wahl, P. and Kern, K. 
Kondo effect of single cobalt adatoms on copper surfaces. Phys. 
Rev. Lett. 88, 096804 (2002). 

19 Nagaoka, K, Jamneala, T, Grobis, M. and Crommie, M. F. Tem- 
perature dependence of a single Kondo impurity. Phys. Rev. Lett. 
88,077205 (2002). 

20 Wahl, P., Diekhoner, L., Schneider, M. A., Vitali, L., Wittich, G. 
and Kern, K. Kondo temperature of magnetic impurities at sur- 



faces. Phys. Rev. Lett. 93, 176603 (2004). 

21 Neel, N., Kroger, J., Limot, L., Palotas, K, Hofer, W. A. and 
Berndt, R. Conductance and Kondo effect of a controlled single 
atom contact. Phys. Rev. Lett. 98, 016801 (2007). 

22 Santoro, G. E. and Giuliani, G. F. Impurity spin susceptibility of 
the Anderson model: a perturbative approach. Phys. Rev. B 44, 
2209(1991). 

23 Pollwein, W., Hohn, T. and Keller, J. Spin-polarization around a 
Kondo impurity. Z. Phys. B 73 , 219 (1988). 

24 Ishii, H. Spin correlations in dilute magnetic alloys. J. Low Temp. 
Phys. 32, 457 (1978). 

Borda, L., Fritz, L., Andrei, N. and Zarand, G. Theory of inelas- 
tic scattering from quantum impurities. Phys. Rev. B 75, 235112 
(2007), 

26 Rubtsov, A. N., Savkin, V. V. and Lichtenstein, A. I. Continuous- 
time quantum Monte Carlo method for fermions. Phys. Rev. B 72, 
035122 (2005). 

27 Wiegmann, P. B. and Tsvelick, A. M. Exact solution of the An- 
derson model: I. J. Phys. C: Solid State Phys. 16, 2281 (1983). 

28 Gubernatis, J. E., Hirsch, J. E. and Scalapino, D. J. Spin and 
charge correlations around an Anderson magnetic impurity. Phys. 
Rev. B 35, 8478 (1987). 

29 S0rensen, E. and Affleck, I. Scaling theory of the Kondo screening 
cloud. Phys. Rev. B 14, 9153 (1996). 

30 Borda, L. Kondo screening clound in a one-dimensional wire: nu- 
merical renormalization group study. Phys. Rev. B 75, 041 307 (R) 
(2007). 



SUPPLEMENTARY INFORMATION 

SI. CONDUCTANCE OPERATOR AND SPECTRAL DENSITY 

Quantum mechanically the fluctuations, or uncertainty, of an observable is related to the variance of its expectation value. 
To obtain an expression for the conductance fluctuations or more specifically the spectral density, one needs an operator for the 
conductance. This can be found by taking the derivative of the linear response expression of the tunneling current, equation (2), 
with respect to the applied voltage, V. Doing so gives 

oo 

G LR M) = ze 2 |T| 2 £ J dt'0(t~t')(t-t'){[C a ,(T,t'),CUv,t)]-[CAr,t),Cl(T,t')]} 

oo 

-ze 2 |T| 2 ^ J dt' &(t-t')(t + t>){[C a (T,t),C a ,(v,t')} + [cUr^),^,^^')}}, (SI.l) 



where C ff (r) is defined in terms of the field operators of the STM and substrate; C CT (r) = ^ , g TM (rcr) , J/ sub (r<7). Although the 
tunneling matrix elements, T, in (SI.l) have different dimensions than that of (1), we will use the same symbol for simplicity. 
Upon taking the expectation value, the second term of (SI.l) generates anomalous Green's functions that only contribute for a 
supercurrent. These terms are responsible for the Josephson effect and will not be considered here. For a non-superconducting 
system the expectation of (SI.l) recovers the well-known expression for the conductance in terms of the local single-particle 
density of states p°(r,uj); (G LR ) Hq = 2vre 2 |T| 2 £ ff pW^KutA, eV). 

The symmetrized spectral density that characterizes the frequency distribution of fluctuations of the conductance about its 
averaged value is defined as 

S(t,u) = - J dte iut ({6G hK (r,t),6G hR (r,0)}) Ho , (SI.2) 
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where SGlr — Glr — (Glr)- This expression for the spectral density assumes time-translational invariance, although 
the full expression for the conductance operator (SI. 1 ) contains non time-translational invariant terms, these vanish for non- 
superconducting systems. Using (SI. 1) and (SI. 2), a lengthy, but straightforward, calculation (SI. 2) gives the zero-frequency 
component of the noise as 

S(t,u = 0) - 87r 2 e 4 |T| 4 [psTM(eF)] 2 X ^ b (r, W = 0) + 327r 2 e 4 |T| 4 [ msTM (e7)] 2 ^ b (r, w = 0). (SI.3) 

Equation (SI.3) is the starting point and motive for the remaining calculations and results of the article. Technically, only the 
vertex corrections of (5a) and (5b) are measured in (SI.3), but in the low-temperature limit, the single-particle contribution to 
these correlation functions is negligible and vanishes at zero temperature. Thus, we will work with the full correlation functions 
for simplicity and clarity. 



SII. GENERATING FUNCTIONAL FOR N-CORRELATED SITES 



In this section we show how one can, in a simple and straightforward manner, obtain arbitrary correlation functions of a 
system consisting of TV-correlated sites in contact with a non-interacting bath. In essence this method reduces all correlations 
to correlations involving only the impurity operators. The remaining impurity problem can then be solved using a variety of 
techniques, such as numerical renormalization group or quantum Monte Carlo. Theoretically the number of impurities N can 
be arbitrary, but in practice it is limited by the complexity of the system at hand, such as the number of bands of the substraight 
and orbitals of the adatoms. In the simplest case, for two-particle properties, an upper bound of four or five Anderson impurities 
could be incorporated using the CT-QMC impurity solver used in this work. 

The partition function Z = Tr e~@ H with respect to the iV-impurity Hamiltonian (6) can be written as a functional integral 
over Grassmann fields (by inserting complete sets of fermionic coherent states) as 

N 

Z =Ii D[b,b]D[c a ,c a ]e- s (SUA) 

a=l ' 

with an action S in Fourier-space 

S = EE E ^LOm) [iu m - (eL - M)] KAiu m ) + E c as (iuj m )(iuj m - E as )c as (iui m ) 

n kcr iu m a,s iuj rn 

+ E E E [ v ^, a sKA i ^ m )c as { i ^ m ) + V£Z as c as (iu m )b%. a (w m )] + s- mt [c as {ibj m ), c as (w m )], (sn.2) 

n,a,s kcr iw m 

where iu> m are fermionic Matsubara frequencies and S- m t is the term of the action (Hamiltonian) that corresponds to the interac- 
tions of the N impurities. The actual form of which is immaterial, as long as it is a local interaction. By introducing Grassmann 
source terms, A and A for each impurity and the bath electrons, the generating functional Za for the system is given by 

ZA= f[f D lb,b]Dlc a ,c a }e~ s expl YY [Aas{iu m )c as {iu) m ) + A as (iuj m )c as (iuj m )] 

ct—1 a,s iuj rn 

+ E E E [KA^ m )bU^J + AU^ m )KA^ m )}\. (SII.3) 

n kcr 

Because we have assumed a non-interacting bath and local interactions of the impurities, the bath electrons remain Gaussian and 
with the identity 



D[r],r]]exp ~JJ dxdx'ri(x)M(x,x / )r/(x / ) + J dxJ(x)rj(x)^(x)+ j dxJ* (x)^(x)rj(x) 
det[M(x,x')]exp\ [[ dxdx' J*(x)^(x)M- 1 (x, x')£(x') J(x') 



(SII.4) 
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can be exactly integrated out of the generating functional (SII.3), leading to 

Z A = det [iui m - (e£ CT - fj,)] f[ J D[c a , cJe~ Soff exp \ EE [A a s{iu m )c as (iu m ) + A as (iuj m )c as (iu m )] 

2^ 2^ 2^ f w _ ( e n^ _ ^) 2^ 2^ 2^ i h] _ _ ,^ 2^ 2^ 2^ 



n,a,s kcr iu;^ 



n,a,s kcr la; 



iuj m - (eL - m) 



where the effective action of the impurities is, in matrix notation, 



(SH.5) 



S c « = -^2^2c^(iu) m )g 1 (iui m )c s ,(iu m ) + S iM [c as (iui m ),c as (i 

s,s f iuJ m 



(sn.6) 



with 



c s (iw m ) 



/ c sa=1 (iui m ) 



\ c sa=N (iuj m ) 



V c sa=N (iu m ) 



c being the transpose of the column vector, and 

[Go 1 ( iuj m )] as a , s , = (iu m - E as )5 a , a >5 s>s > - A aa>a > s >(iuj m ), 
where the hybridization or Weiss field is 

A QS ,<i' S '( iw m) = E 



n.kcr 



V£ n , , 

kcr, as kcr, a s 



The determinant that appears in (SII.5) can be analytically evaluated, but it is actually never referenced in a calculation of a 
correlation function, as it is canceled by an identical term appearing in an overall normalization factor. 

Using (SII.5), arbitrary correlation functions can be obtained by suitable functional differention. As can be seen, bath cor- 
relators (apart from trivial non-interacting terms) are found from impurity correlators simply by attaching "tails" of the form 

iui _^ e '° s _ A1 ) or iuj _( e n s _^ ■ For example the single-particle Green's function of the bath is 



Gn,n / ■ ■ \ 



1 



A=0 



/"."'-EE; 



V, 



kcr, as 



a., a' s.s' 



( C a,s(^m)<V, S '( iw m')) 



k'cr' ,<x' s' 



(sn.7) 



which is just the multi-impurity generalization of the well-known T-matrix expression, commonly found by equation of motion 
methods. Similarly the two-particle Green's function of the bath is given by 



kieri,k 2 <T2,k3Cr3,k4er4 



/~r"l,«2 ( 
kieri,k 2 <T 2 



(1^,^2,1^3,^4) 

n 3 ,n 4 



// ■ .!)■? / - \ 



E 



K" 1 

kicri,aisi 



k 2 a 2 , Q2S2 A". 1 '" 2 !?;" 4 (iw l5 kj 2 , iu> 3 , k3 " 3 ' a3S 



y*„ 4 

k,j (7.| .( t 1 .S ; 



S'2 ,^3,54 



(SII.8) 



where 



E 



fll ,«i ( lUJ l) G s 2 ,s' 2 (^2)r s /^ 4 3 s / 4 («^1, *W 2 , IW 3 , *W 4 )G S , 3 S3 («^ 3 )G S ,% 4 (tWj 



1' 2' 3' 4 
s l ' s 2 ' S 3 ' S 4 



(SII.9) 
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is the non-trival part of the impurity two-particle Green's function, i.e., 

+ E (S11.10) 



here T is the reducible vertex of the impurity problem. 



